Dark Energy Content of Nonlinear Electrodynamics 
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Quasi-constant external fields in nonlinear electromagnetism generate a contribution to the 
energy-momentum tensor with the form of dark energy. To provide a thorough understanding 
of the origin and strength of the effects, we undertake a complete theoretical and numerical study of 
the energy-momentum tensor T^" for nonlinear electromagnetism. The Euler-Heisenberg nonlinear- 
ity due to quantum fluctuations of spinor and scalar matter fields is considered and contrasted with 
the properties of classical nonlinear Born-Infeld electromagnetism. We also address modifications 
of charged particle kinematics by strong background fields. 
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I. INTRODUCTION 

Recent analysis has constrained the dark energy to 
have almost exactly the characteristics of a cosmologi- 
cal constant: equation of state w =_p/p « — 1 and spa- 
tially homogeneous distribution P, Q . This means that 
the dark energy is present in the energy-momentum ten- 
sor proportional to g^'^ as has often been discussed in 
the context of vacuum energy Q . Proportionality to 5^'' 
docs not exclude that dark energy originates in proper- 
tics of ponderable fields and matter; a nonvanishing trace 
of the energy-momentum tensor Tj^ = g^'^T^^, where 



T — 



(1) 



signals presence of a dark energy term in T^^. 

A good candidate for dark energy is a physical the- 
ory that, without interactions, has a tracclcss energy- 
momentum tensor. Extensions of this theory including 
interactions can then be analyzed for the generation of 
an energy-momentum trace, and by comparison with the 
bare theory, the physics of the energy-momentum trace is 
more easily extracted. An effort in this direction has been 
made in the context of vacuum structure of quantum 
chromodynamics (3], we show here that standard quan- 
tum electrodynamics (QED) has similar features which 
are more easily accessible. 

Our point of departure is thus Maxwell electromag- 
netism, whose energy momentum tensor. 



^Max 



(2) 



is explicitly traceless. Even in absence of external sources 
the Maxwell field equations are incomplete due to in- 
teraction with the electron-positron vacuum fluctuations 
which are present at the length scale Ae = H/rrieC. 
At distances of comparable magnitude (A ~ Ae) these 
arc vacuum polarization effects which impact precision 
atomic physics experiments. For long distance (A 3> Ae) 
one obtains the nonlinear effective theory of the photon 



studied in depth by Euler, Kockel and Heisenberg, and 
Schwinger 

This nonlinear Euler-Heisenberg (EH) theory of elec- 
tromagnetism is just one of many possible effective ac- 
tions. Beyond the EH-QED framework, we can imagine 
writing down a more complete theory containing all ef- 
fective interactions, reducing in the long wavelength limit 
at the classical level to Maxwell's equations. Born-Infeld 
(BI) theory, designed to regulate point-particle-induced 
divergences , can be thought of as an effort to provide 
a more complete theory of electromagnetism, though it 
too is now considered as an effective theory arising 
from string theory. 

We show how in nonlinear electromagnetism, the clas- 
sical energy-momentum tensor Eq. ^ is modified by two 
quantities: a dielectric function e, which scales the con- 
tribution of the Maxwell energy-momentum in the total, 
and the trace of energy momentum tensor T^, which is 
a new contribution. We have proposed that arising 
in QED can in this regard be viewed as the modification 
of the vacuum energy by the presence of electromagnetic 
fields 3 , and indeed the connection of dark energy with 
a form of vacuum energy has been discussed before 0, Q ■ 
In this work we fully develop this conjecture, evaluate the 
dielectric function and compare the behavior of EH and 
BI nonlinear theories. 

As is well known, a theory with a traceless energy- 
momentum tensor is invariant under scale changes. 
Maxwell's electromagnetism is the prime example: the 
energy-momentum tensor is traceless and indeed the clas- 
sical theory of radiation is scale invariant. In a gauge 
theory, the related conformal symmetry can be sponta- 
neously broken, the study of which has a long and dis- 
tinguished history, originating with the rise of quantum 
chromodynamics (QCD) as the theory of strong interac- 
tions [id . [Tlj . We will rederive here some results of the 
extensions of these studies to QED [l^, [13 . 

Scale invariance can be broken by the explicit appear- 
ance of a dimensionful quantity in the theory, such as 
the mass of the electron in QED. In fact, any nonlinear 
electromagnetism requires a scale with the dimension of 



2 



an electrical field cEq, in order to render the Lagrangian 
dimensionally consistent. We express this scale in terms 
of a mass M: 

sEq = — 

he 

For the Eulcr-Hcisenbcrg (EH) effective action the intro- 
duction of M is a natural step to take as the nonlinearity 
is of quantum origin and M ~ m/y^, where a = 1/137 
is the usual fine structure constant. For the BI theory it 
is a matter of convenience to use mass rather than length 
as the scale, converting one into another using ft. 

If indeed BI theory is a weak-field limit of string theory 
in which the nonlinearity is a consequence of high mass 
quantum fluctuations [8| the appearance of ft would be 
appropriate, and the associated scale could be as large 
as the Planck mass, Mpi = 1.2 x 10^^ GcV. It should 
be noted that current experimental limits as well as EH 
nonlinearity probes a scale below 100 MeV, thus a string 
related BI theory maybe quite removed from the present 
experimental reality. We compare EH and BI theories 
mainly because their behavior is very different. 

In QED and BI, the presence of scale which breaks 
the conformal symmetry is explicit, and the energy- 
momentum trace is not "anomalous," unlike the case of 
QCD. The scale of the nonlinearity is, as we shall show, 
the determining factor of T^. This fact is the simple yet 
important and original theoretical observation presented 
in this paper. Having thus suggested the inter-connection 
of dark energy, conformal symmetry and the presence of 
scale in the theory, we leave issues specific to conformal 
symmetry to future work and here focus on the physics 
of and its origins in nonlinearity of the theory. 

We derive in Section|TT]the field energy-momentum ten- 
sor and explicitly connect its trace to nonlinearity of the 
electromagnetic theory. To compare relative magnitudes 
and suggest new constraints on a Born-Infeld-type com- 
pletion of electromagnetism, in Section lllll we evaluate 
the BI modifications to the Maxwell energy-momentum 
tensor. In clarification of conflicting statements present 
in the literature, we begin our discussion of quantum 
electrodynamics in Section IIVI with a new derivation of 
the relationship between the electron-positron conden- 
sate and the energy-momentum trace based on the ex- 
plicit origin of the trace in nonlinearity of the theory. Ex- 
tending a technique of resummation of the action intro- 
duced by Miiller, et al. [l^l, we then provide complete nu- 
merical evaluations of the condensate, energy-momentum 
trace and dielectric function for QED and spin-0 quan- 
tum electrodynamics. These evaluations display strik- 
ing analytical features not before apparent in the Eulcr- 
Heisenberg functions. We compare BI and QED contri- 
butions to and show that QED vacuum fluctuations 
remain dominant given the experimental constraints. 

In the final Section |V] of this report, we discuss the 
kinematics of charged particles moving in external fields. 
The vacuum of a nonlinear theory is studied as a pon- 
derous medium with nonlinear response. The Lorentz 



force is preserved, but the breaking of the superposition 
principle results in effective potentials for charged parti- 
cles moving in external fields that are not automatically 
obtained from the Lorentz force. 



II. ENERGY-MOMENTUM TENSOR OF 
NONLINEAR ELECTROMAGNETISM 

Setting from now on ft = c = 1, we can consider the 
effects of Eq ■= M^/e mass Af or length I ~ 1/M scale, 
where M can be as large as a string theory scale or as 
small as the mass of the electron. The consequences are 
best seen writing the effective action in the form 

^---'^ + ^-^'^(i|i']^).^o.-'^ (3) 

presented here clS 0j function of the Lorentz scalar and 
pseudoscalar 

S (4a) 

V := \f:^F-^=E-B. (4b) 

As noted, classical, linear electromagnetism must consti- 
tute the limit of Eq. fS]) for fields small as measured in 
units of Eq, and only the classical theory docs not require 
a dimensioned scale. 



A. Dielectric function and Trace 

To understand the implications of the dimensioned 
scale we consider the explicit form of the energy momen- 
tum tensor ([1]), separating the tracelcss Maxwell part. 
For a general function Voff (5, 7^), we obtain 

= (-^)(,9M.5-f^^AF/) (5) 

Comparison with Eq. ([2]) shows the energy-momentum 
tensor of Maxwell's electromagnetism is modified by a 
dielectric function, —dVc^/dS, to be discussed below. 
Using Eq. ([3]) we simplify the second term to 

This form is very useful because it provides a simple 
means of calculating the trace directly from the effective 
action. The importance of Eq. ^ lies in its distillation 
of the physical source of conformal symmetry breaking 
in any nonlinear theory of electromagnetism: 

Terms linear in the invariant S cannot contribute to 
the right side of Eq. © since they cancel explicitly on 
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the left side of Eq. ([6]). Such contributions must therefore 
be omitted from Vgs in the study of energy-momentum 
trace, and hence we have introduced the barred fos to 
denote the nonhnear components of the effective poten- 
tial. Letting V^P denote the remaining linear terms in 
the Lagrangian, we have the decomposition 



(7) 



The lowest power in 5, 7^ is 2 due to preservation of both 
parity and charge conjugation symmetry, which respec- 
tively require an action even under parity transforma- 
tions and even in the coupling e, hence even in the field 
strengths (already true of any Lorentz scalar). A nonzero 
imaginary part of the action entails breaking of time re- 
versal symmetry, though the trace Twill be small under 
laboratory conditions. These symmetry arguments imply 
that for field strengths below the critical scale E ^ m^, 
the energy-momentum trace must be at least 4th order 
in the fields. 

Solving the partial differential equation 



cS 



S 



dS 



(8) 



displays one obvious class of nonlinear Lagrangians that 
have traceless energy-momentum tensors, namely 
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(9) 



Comparison with Eq. ^ reveals the reason: such La- 
grangians are conformal. Since this class is non- 
perturbative in at least one of the field invariants and 
we are interested in having an energy-momentum trace, 
we exclude theories of the form in Eq. ^ as well any 
other non-perturbative actions that may satisfy Eq. ([8]), 
despite their inherent interest. 

In view of Eqs. (j5l6|) . we summarize 



the quantum effects and their connection to nonlinearity 
of the action. 

The alternative and equivalent representation, often 
used in the electromagnetism of nonlinear media (for ex- 
ample see section 8 of [l6j). 



is using the displacement tensor 



dC 



dFn. dS 



dV 



(11) 



(12) 



The trace is now distributed into several components 
r = W^F^^ -AC = E-D + B-H -AC. (13) 

The origins and properties of T are obscured by the 
constitutive relation H^'^ {F^^^ ,F^^). Our expression 
Eq. (fTO|) evidences the departure from the classical theory 
more clearly in the context we consider. 

B. Stress-Energy Density 

Some further notable properties of the stress-energy 
density of the nonlinear EM field will be collected here. 
The energy density in the frame of the metric, i.e. the 
quantity entering Einstein's equations, is 



^00 



1 



-T. 



The EM-stresses T'^ have the same structure as in the 
Maxwell theory, as is already evident from the format of 
Eq. dUD: 



Ma 



5in -T 



(14) 



e = —■ 



dS ' 



Tt^ = T = 



(10a) 
(10b) 



Tf^"" the Maxwell 



Here e is the dielectric function 

energy-momentum tensor, and T the energy-momentum 
trace. Interestingly, T/4 ^ A/2 provides a dark energy 
or Einstein-like cosmological constant, while the trace- 
less part is the same as in Maxwell theory, up to the 
multiplicative dielectric function. 

Eq. (jlObp though derived here for the case of electro- 
magnetism has a much wider domain of validity. The 
implications of the separation in Eq. (|10p have been pre- 
viously noted in context of the photon propagation ef- 
fects [13: but Twas not given in the form Eq. (jlObp . 
nor has in any form the trace T been computed. We 
shall below demonstrate how the concomitant identities 
Eq. (jlObp provide new physical insight in understanding 



The trace T acts to compensate the forces T^^x tear- 
ing the field sources apart in Maxwell electromagnetism. 
For this reason, for example in BI theory the T'^ van- 
ishes allowing a stable charged particle without material 
stresses. A sufficient condition for this to be true is that 
the point particle solution satisfies lim^^o r^T^'^ = [l3]. 

Because the energy momentum tensor is conserved we 
have T^y^^j = dT^^/dx'^ = 0, which is a covariant relation 
true in any frame. A differential conservation law leads 
to an integral conservation law by integration over the 
observer's hypersurface: 

J = 0, or J d^ayTf"" = J d^ayT^"', (15) 

i.e. the energy-momentum flow through surface 1 is the 
same as later through surface 2. It is common to choose 
an observer at rest in laboratory so that d'^o^ — u^d^x, 
with u"^ = (1,0,0,0). 



4 



For nearly homogeneous fields we can omit the 3- 
volume and consider a conserved 4-monientum density 
of the EM field 



fL rpl-LV 



finding the well known result for the rest energy and 
Poynting vector of the classical field. This result is easily 
generalized to the nonlinear electromagnetism: 



e ^ + _,eSxB 



(16) 



showing the appropriateness of calling e the dielectric 
function, since it plays the role of the dielectric constant 
E = eD when considering electric charge in vacuum. For 
Maxwell electromagnetism e = 1 and T = 0. 

It is an elementary exercise to show that in the Maxwell 
limit the proper energy density, or its "mass density," is 



„Max u 



{E^+B 
4 

(£2 - B^ 



2\2 



{E X Bf (17a) 
{E ■ Bf=S^ +r^il7h) 



Generalizing to nonlinear EM theory, we find the local 
mass density of the field to be 



T/4 ^ A/2 provides thus both a 'dark energy' and a 
'mass density' of the electromagnetic field. 



where the particular combination of S, P terms derives 
from the extension of the space-time metric with the an- 
tisymmetric field tensor, Fp_u. In the weak field (infinite 
mass) limit Maxwell's theory indeed arises, 

y(BI) ^ _S. 

For the Born-Infeld case the dielectric function is 



(BI) 



^BI 



cff 



dS 



[l + 2S/M^ ~{V/M^f]-^''^, (22) 



which exhibits a formal analogy to the 7-factor familiar 
from special relativity, though with two different limits 
as 5 or 7^ respectively approach the limiting value 
(sec figure [TJ left). The dielectric function goes over from 
suppression (e < 1) to augmentation (e > 1), when the 
magnetic component of the field becomes subdominant, 
which corresponds to crossing the line 25 = from the 
lower right. 

As presented, Eq. p9p is in the form required by 
Eq. ([3]). However, the BI action contains no terms lin- 
ear in S; T is identically —M{dVcs/dM). We obtain for 
the BI energy-momentum trace using the relation Eq. ([6]) 



BI) 



4M^{eBi{l+S/AlH) - 1), 



(23) 



4M^ 




which in the latter form is manifestly positive-definite, 
just like the cosmological constant. For small fields we 
expand the second form in Eq. (|23p to obtain 



III. BORN-INFELD ELECTROMAGNETISM 

As demonstrated in the preceding discussion, an in- 
trinsically nonlinear theory of electromagnetism (in most 
cases) entails an energy-momentum trace, and we begin 
by studying a non-quantum example of a nonlinear al- 
ternative to Maxwellian theory. Historically, Born-Infeld 
electromagnetism was introduced in order to solve the 
infinite self-energy (and self-stress) problem of a point- 
like electron arising in consideration of the radial electric 
field Er of a point charge q: 



U = I £x-El 



00, for Er = 



To remedy this, Born and Infeld took inspiration from 
special relativity, considering the action (note that we 
follow the modern convention, opposite in sign to the 
original paper Q, and also, recall remarks about the scale 
M above Eq.©): 



(BI) 



off 



M'^ (1 - v/1 + 2S/M^ - CP/A/4)2) (19) 
^-\/^), (20) 



M 



h = det h 



5.. + ^, (21) 



7-(Bi) 



M4 ^i + 25/A/4 



(24) 



In figure [T] we show the dielectric function and energy- 
momentum trace for strengths up to the maximum field 
strength. The fimctional behavior of both is smooth, 
though & S,V functional asymmetry develops at large 
values of the fields. 

In contrast to the classical theory analyzed here, we 
note a recent report suggesting that the quantized BI 
theory as studied on the lattice may be conformally sym- 
metric list. 



IV. EULER-HEISENBERG 
ELECTROMAGNETISM 

The Euler-Heisenberg effective action is well known: 



for Dirac fermions, and 



(1 — eas cot(eas) e6s coth(e6s)) (25) 



* eft 



0+5 



IGtt^s^ 



(eas csc(eas) e6s csch (e6s) — 1) (26) 
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for charged scalars, in which 

a \/ y/S^ + andfe \/ y/S^ + V^- +S. (27) 

The characteristic strength of fluctuations in the matter 
field is made explicit in the appearance of m which is the 
mass of the matter particle which has been integrated 
out, in QED it is the mass of the elec tron. 

For vanishing V we have a — S and b 

y/\S\+S. Thus when £: wc find a ^ 0, 6 ^ \B\ and 
when B = we find a ^ \E\,b ^ Q. In this sense, b plays 
the role of the generalized (Lorcntz invariant) magnetic 
field and a that of the generalized electric field. The con- 
stant subtraction ±1 removes the (divergent) zero point 
energy of free electrons and positrons. The difference 
in normalization reflects the doubling of the number of 
degrees of freedom for spin- 1/2 particles, and the over- 
all sign corresponds to the difference in sign of vacuum 
fluctuations between bosons and fcrmions. 

A. The condensate {^pip) and trace T 

Equation (jlOb[) provides a direct means of calculating 
the energy-momentum trace, but its connection to vac- 
uum structure in QED and its deformation by the ap- 
plied fields which induce the EH nonlinearity is encoded 
in m{dVcff/dm). Consider the Feynman boundary con- 
dition Green's function of the fluctuating matter field in 
presence of the electromagnetic field, which determines 
the vacuum fluctuations, 

^91offO^^ — ) _ [Spix+e, x—e, to) — S'p (28) 

dm e-fO 

Here e is a time-like vector, and S'p^ is the free field 
Feynman Green's function. In general: 

SF{x,x') = -t{T{ij{x'Wx))). (29) 



We rewrite the right side of Eq. ([^5]) using the elemen- 
tary form of Wick's decomposition theorem 

T(^(x')^(x)) =: ^{x'Mx) : +(0|r(^(x')^(x))|0), 

where the normal ordering is with respect to the 'no-field' 
vacuum. Taking the expectation value of this relation at 
a single space-time point in the 'with-field' vacuum |) we 
find 

(: ib^ :) = {\Ti^ix)^{x))\) - (0|r(7A(x)^(x))|0) (30) 
= iSpix, x) — iSp{x, x) 

with the same e-limit as in Eq. (|28p implied for equal 
propagator arguments. Usually the trace is implied by 
commuting the fields and the normal ordering symbols 
omitted, 

to(: Tpip :) —m{i]jij}), 

hiding the important operational definition Eq. (j30p of 
what we now recognize as the fcrmionic condensate. 
Thus, we see that the condensate derives from the dif- 
ference of normal ordering in the no-field (also called 
perturbative) vacuum and the with-field vacuum. Fur- 
thermore, one must keep in mind that the subtraction of 
the unperturbed vacuum term follows directly from ap- 
plication of the rules of QED and is not a consequence of 
arbitrary removal of zero-point energy. For this reason 
our discussion has no bearing on the zero-point energy of 
quantum field theory and/or its gravitational coupling. 

Equations (HH), and then combine with the 
result 

TO— = — to(-0'!/'), (31) 

dm 

a known and widely used relation in non-perturbative 
QCD. We will evaluate it for the case of electrons in pres- 
ence of external electromagnetic fields. 
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We undertook the derivation of the Fermi condensate 
in the preceding section in order to emphasize that the 
condensate is not in general the energy-momentum trace, 
though there is a relationship. Using Eqs ([7]), (jlOb|) . and 
(|3T]) we find 



r 



If and only if V^^^ is just the action of classical electro- 
magnetism, the first term on the right hand side, linear 
in S, vanishes. 

However, the coupling to a fluctuating quantum field 
complicates the issue significantly in that it produces a 
contribution linear in S with a coefficient which is a func- 
tion of m and renormalization scale. It is the logarithmic 
divergence of Eqs. and which in the limit (5 — > 

is also proportional to S and thus appears in V^^\ The 
standard procedure is to absorb the divergence into the 
definition of the charge in the process of charge renor- 
malization. Eqs. (|25|) and ((26)) assume a cutoff regulator 
which of course also introduces a scale. Alternatively one 
can use dimensional regularization, the spinor effective 
potential can be written 



(1) 



dm 



(32) 



The essential relation Eq. ([36)1 must be preserved at 
any number of loops in the effective action in its suitable 
generalization. In particular, if the condensate were eval- 
uated to two loops, the coefficient of the first term must 
become the two-loop /3-function, i.e. 



2a 
3^ 



/?(«) 



2a 
3^ 



a 

2^ 



+ 



(37) 



Our result Eq. is not obvious if one evaluates the 
energy-momentum trace after renormalization has been 
carried out. The logarithmically divergent term in V^s is 



(1) 



62- 



eff 



87r2 



ds 
— f 

0-1-5 S 



—SMm^S) 
in 



(38) 

in which 5 = f/M^, some large mass or momentum ex- 
traneous to QED. Before presenting the EH action, this 
term is absorbed in the process of charge renormaliza- 
tion. To restore its contribution one must realize the 
mass dependence of charge renormalization. The relation 
of Eq. (pSj) to the QED (3 function [Tl| demonstrates why 
use of the renormalized Vch with the incorrect identifica- 
tion T = mdVcfi I dm leads to the correct result as shown 
in Eq. dMl), just as it was developed for QCD [i,[l3,[lH- 



* en 



^ds e 



2 „3-c 



(1 — eas cot(eas) e6s coth(e6s)) (33) 



and similar for the charged scalar case. 

For any finite e, Eq. ((33)) is finite and can be dif- 
ferentiated with respect to m, and hence the quan- 
tity m{dV^g / dm) is finite, allowing the considera- 
tion of a vanishing e. In particular, the differenti- 
ation midVcs/dm) renders the lowest order contribu- 
tion m{dV^^ /dm) finite. The condensate and energy- 
momentum trace are thereby independent of renormal- 
ization procedure (as they should be), and obey the non- 
trivial relationship expressed by Eq. ([5^ . The finiteness 
of T for fermions is also discussed at length by Adler et 
al [13. 

Separation of the term linear in S in the Fermi case 
shows 



(1) 



off 



e2 62 



dm 



87r2 



ds 



2a 
3^ 



S. (34) 



Thus for the EH effective action, the decomposition in 
Eq. ([7]) becomes 



m , = m 
dm 



dm 



(35) 



and in turn with emphasis on properties of the vacuum: 



2(T 

r' = — (5)+m(^V). 
3tt 



(36) 



Equation (|5S)) corresponds to Eq.(2.17) in [T3|. It is of 
importance to note that there is considerable cancellation 
between the two terms on the right hand side Q. 



B. Properties of T in QED 

We obtain the explicit form of the condensate in ex- 
ternal fields combining Eq. (|3T)) and Eq. ([25)) 



m{ip^p) = m" 



(eascot(eas) x (39) 
ebs coth(e6s) — I). 



The condensate vanishes in the absence of field, as it 
should, since a; cot a; —i- 1 and xcothx —f 1 for x — > 
0. The term quadratic in the fields has been discussed 
above, it must be subtracted to arrive at the integral 
representation of the energy-momentum trace: 



47r2 



{eas cot(eas) x 



(40) 



e6scoth(e6s) - f (fe2 - a2)s2 ) . 

3 



To study the integrals Eq. ([55)) and Eq. ([HI)), we intro- 
duce a transformation that will be helpful in dealing with 
the non-analyticities generated by the electric field. The 
detailed calculations are carried out in the Appendix [A) 
and we simply state here their results. The integral rep- 
resentations obtained along the way have better conver- 
gence properties than Eq. ([39)) or Eq. ([40)) , particularly at 
strong fields B,E ^ Eq. 

In a magnetic background field, the condensate can be 
written 



m{'ip'il}) 



27r2/3' 



ln(f 



-f3's 



ds, 



(41) 
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where /3' := nrn'^/eB = ^/{B/Eq), with Eq = m'^/e. 
The energy-momentum trace is obtained by removal of 
the leading term quadratic in the field (recall Eq. (j36|) ). 
with the result that 



'ln(l - e-'3'^) 

rT^2 ' 



(42) 



is manifestly positive definite. 

For the electric field, the poles are resummed into a 
logarithmic winding point, with the result 



27r2/3 



ln(f - e-^'^) 



1 



ds 



(43) 



where /? := tt/{E/Eq). The trace in the pure electric 
background is 



°°s^ ln(l-e-^« 



ds. 



(44) 



Equations ([321) and ([331) could cquivalently be obtained 
by taking B iE in the respective magnetic expressions. 
The condensate behaves as (ei?)^ for small fields, but 
the poles displayed in Eq. (|A7p give the condensate and 
the energy- momentum trace a nonzero imaginary part, 
which can be evaluated from Eq. ([321) or ([331) recalling 
the identity 



PV- 



't7ri5(a;). 



As befits its role contributing to the proper mass of the 
nonlinear electromagnetic field (recall Eq. ([TSl) ). 



ImT = — m — ^ > — e 
47^2 ^ n 

n=l 



47r/3 



ln(f 



(45) 



is manifestly negative and strongly suppressed for field 
strengths less than 0.1£^o- This is consistent with direct 
differentiation of the positive imaginary part of the action 
V^g evaluated by Schwinger 



{eEf 



Stt^ ^ — ' n 

n—l 



oo 

Yi2 



Observe that m{ip^) in external magnetic fields is 
negative, as is manifest in Eq. (|41[) . while the energy- 
momentum trace, Eq. (j42|) is positive. With the perspec- 
tive that the energy-momentum trace represents the en- 
ergy arising from deformation of the vacuum, a negative 
value would imply that the vacuum state is unstable, for 
example to the spontaneous generation of strong mag- 
netic fields if T were identically m{ip'ip). 

The energy-momentum trace in general combinations 
of electric and magnetic can also be cast Eq. (jA9[) remi- 
niscent of our prior study of the special cases of electric 



and magnetic fields. Using b in the definition of /3' and a 
in /3, the numerically useful representation is: 



(46a) 



k=l 



ds ! 



1 



E^^coth^. (46c) 
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A further rcsummation (see appendix [^ displays a more 
statistical form of the integrands: 



la = -^E 

h = aE 



>o s i In f if^l 

e/'^" - 1 
s — Arctan s 



e 



1 



■ds 



ds (47a) 
(47b) 



where a = ±f and n = 0, f,2... in accordance with the 
Landau levels apparent in the quasi-Hamiltonians 

m'^na = rn^s + (2n + f - a)eh (48a) 
m^li!„ = rn^s + (2n + 1 - a)ea. (48b) 

The slower convergence of la (h) for b/a <^ I {a/b 1) 
means Eq. ([47]) is not as advantageous as Eq. ([46l) for 
numerical work. These expressions clarify the electric 
and magnetic field limits and thereby the correspon- 
dence between figures[2]and[3]below when one also recalls 
z coth z ^ 1 for z ^ 0. 



C. Numerical Evaluation of QED T 

We evaluate numerically the condensate and the 
energy-momentum trace for arbitrary constant, homo- 
geneous background electromagnetic field. We consider 
first a magnetic only background field in Eq. p9|) , a form 
which is easily integrated numerically. The behavior is 
displayed in figure [D as the solid upper (red) line. Not- 
ing that f{x) = (xcothx — l)/x'^ is essentially con- 
stant for a: < 1, we recover the quadratic dominance 
~m{'tpip) oc [gBY for small fields. Figure [H confirms the 
small field (mass-dominant) quadratic behavior and large 
field (field-dominant) linear behavior The results 

of numerically integrating (|43p also appear in figure 
but since the result is opposite in sign to the magnetic 
field case, we show the negative of the result as the lower 
(black) curve. 

Eqs. ([331) and ([311) are plotted on the right in figure [H 
and interestingly, wc sec that at _B = 9£'o T changes sign 
from positive at low fields to negative at high fields. This 
result is not apparent in the perturbative expansion, but 
can be understood from the Cauchy-Riemann equations 
in view of the rapidly changing imaginary part. Hav- 
ing the appearance of a tt phase flip, the feature may be 
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FIG. 2: The trace T"^ normalized by the electron mass for spin- 1/2 matter fields in magnetic (left) and electric (right) fields, 
evaluated perturbatively (dash-dot line) from the first two terms in Eq. (|50|) and exactly (solid line) from Eqs. (|42|) and (|44[). 
At left, the condensate in magnetic-only backgrounds is included for comparison of magnitudes, and at right, the imaginary 
part present for electric fields is included. The trace changes sign from positive to negative at E = 9Eo and for higher fields 
the negative {—Tj^) is plotted. 



related to the rapid dissolution (attosecond timescale) 
of fields at magnitudes surpassing Eq [2^. The meth- 
ods developed here are not suited to the dynamics im- 
plied by such fields and the processes used to obtain 
such field strengths. For comparison, we also plot the 
weak field expansions, derived from the original EH ex- 
pression with the power-law semi- convergent expansions 
of cot and coth. The first two terms for the condensate 
and energy-momentum trace are 

,7,1 m'^ E"^ £^ , 

_^(^^) (49) 



«10"' 



127r2 El 



907r2 E^ 



(50) 



907r2 E^ 3157r2 El 

where f ^ = 52 £2 _ _j^2 magnetic or electric 
fields, respectively. 

The condensate and energy-momentum trace for more 
general field configurations are evaluated using the 
rapidly convergent sums in Eq. (|A9|) . and the results 
displayed in figure [3] An examination of the conden- 
sate for arbitrary field configurations, which does include 
the term linear in S (but is not displayed here), shows 
nontrivial features near the line S ~ V. This is re- 
fiecting on the expectation [l^, of the infiuence of zero 
modes, known to be present when the field is self-dual, 
i.e. S ^V. Comparison of the appearance of figure [3] 
with the BI result, figure [T] reveals a profound difference 
in these results, with the vacuum fluctuation effective ac- 
tion being more 'edgy', and suggesting the possibility of 
a more singular behavior in the full multiloop strong field 
case. 

The weak field expansions in general electromagnetic 
backgrounds give 
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FIG. 3: The spinor energy-momentum trace for general E,B 
fields, parameterized by the Lorentz invariants, S,V as de- 
fined below Eq. (|3]). As only appears in VcS, we plot only 
positive V. Although the figure has a logarithmic scale, the 
trace crosses zero at the dark (grey) line, going from positive 
to negative for very large electric fields <S <S — 1, as seen also 
in fig. [2I (Color online.) 
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(52) 



using S,V normalized to the natural field strength, Eq. 
We can compare perturbative EH result ([52]) with the 
Born-Infeld weakfield T Eq. ((23|) , but because the coeffi- 
cients of 52 and 7^2 jo not match in QED, we obtain two 
values for the Bl-type limiting mass, in obvious notation: 



Ms 
M-p 



*/#7TO = 15.16m, or 



(53) 



4/ 45 

28ce 



m = 13.18m. 
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D. T of Scalar QED 

The Eulcr-Hciscnbcrg-Schwinger calculation of the ef- 
fective action is easily extended to the case of a charged 
scalar field. The spin-0 particle has no Pauli spin cou- 
pling UfiuF^'^ in the Hamiltonian, and the proper time 
integral is evaluated for the Klein-Gordon equation with 
covariant derivative, 



[D'^ + m^)(j) = 0, 



2A„ 



(54) 



leading to Eq. (|26|) displayed above. The condensate 
takes the slightly different form, (</)(/>*). Analogous 
manipulation of the proper time representations as above 
(see section lTV A[) leads to an identity similar to (|3ip : 



dm 



In turn, we have 



-s 



(55) 



(56) 



which displays for the scalar field the large cancella- 
tion between the photon and matter condensates. Then, 
since the explicit form of the energy-momentum tensor, 
Eq. ([5]), remains valid, we need only manipulate the ac- 
tion of scalar electrodynamics. 

The alternating sign in the mcromorphic expansions 
of the CSC and csch functions (see appendix \^ leads 
upon partial integration to the opposite sign, or opposite 
"statistics" in the logarithm ln(l ± e~'^*), as was already 
remarked upon by |2Q] . Thus, we have for magnetic only 
background fields 



47r2/3' 



ln(l 



,-13's 



ds 



and for electric only fields 

4 

m 



ln(l 



-0s\ 



-ds 



(57) 



(58) 



47r2/3 7o 1 - s2 

Results of numerically integrating equations (j57p 
and ([55]) appear in figure [H As before, the condensate 
appears with opposite signs when comparing the electric 
and magnetic backgrounds, and the pole structure 
remains consistent for electric and magnetic fields, irre- 
spective of particle type. The imaginary part induced 
by the electric background similarly becomes 



1mm- 



dV: 



off 



dm 



, eE 



8TTf3 



ln(l+e-'3), 
(59) 

which is again negative in continued agreement with the 
role of T in the proper mass-energy of the nonlinear elec- 
tromagnetic field. 

Finally, the energy-momentum trace generated by the 
scalar field quantum fluctuations is 



47r2/3' 



■ ln(l 



-I3's 



ds. 



(60) 



in the magnetic background, and 



' ln(l + 6-'=*'') 



ds 



(61) 



in the electric background. Similarly positive for all but 
the highest electric fields, the scalar energy-momentum 
trace is exhibited in figure S] on right. The effect of the 
fermionic statistics is apparent in the shift of the zero 
crossing to E ~ 2Q.2Eq. In the figure, we also compare 
the weak field expansions. 
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(62) 



(63) 



where £^ = or £^ = —E"^ for magnetic or electric 
fields, respectively. 

The results for general electric and magnetic back- 
grounds are 



ds s 



1 — s2 — ie ^ 

k=l 



^ csch 



fc/3 



h = 



'dss^ -^(^ife-'^P'Trkp' nk(3' 

csch ■ 



1 



52 

k=l 



k(3' 



(64a) 
(64b) 

(64c) 



which is plotted in figure [5) The statistical form anal- 
ogous to Eq. P7|) can be obtained from that expression 
by changing to the fermionic -I- sign in the denomina- 
tor and omitting the spin sum and a term in the quasi- 
Hamiltonians Eq. ((15]) : the spectral fimctions 



2 



l-x 



1 



and s — Arctan ; 



remain unchanged. The weak field expansions for general 
electromagnetic backgrounds are 



3607r' 



247r2 
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6307r2 



(65) 



(66) 



using S, V normalized to the natural field strength, 
Eq. Again, we compare with the Born-Infcld weak-field 
energy-momentum trace (see Eq. (|23p ). and as above, the 
coefficients of 52 and are not the same. For the scalar 
quantum theory then, the two corresponding values for 
the Bl-typc limiting mass are 



Ms 
M-p 




18.64m, and 
30.32m. 



(67) 
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FIG. 4: For spinless particle fluctuations, we display the the energy-momentum trace in units of for both magnetic-only 
(left) and electric-only (right) backgrounds. Again, the condensate is displayed for comparison in the magnetic case, and the 
imaginary part of the trace present in the electric background. Note the change in sign of the trace is now ai E — 5.5Eo. 
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FIG. 5: The scalar energy-momentum trace for general E,B 
fields, Eq. (|64|l . parameterized by the Lorentz invariants. The 
dark (grey) vertical line indicates the change in sign of the 
trace to negative values for very large electric fields. The 
transition is present up to arbitrary values of V as in fig. [3] 
(Color online.) 

E. Euler-Heisenberg Dielectric Function 

As we will discuss below, the kinematical and gravita- 
tional effects of a trace contribution to the electromag- 



netic energy-momentum differ strikingly from the clas- 
sical Maxwell energy-momentum Eq. Although this 
fact should make experimental verification of the pres- 
ence of the trace term easy in principle, the relative 
strength of the classical contribution cannot be ignored 
in the study of real physical systems. Thus, wc complete 
the analysis of the energy-momentum tensor of Euler- 
Heisenberg electromagnetism with evaluation of the di- 
electric function e = —dVcff/dS. 

The EH actions are corrections to the classical —S, so 
the total dielectric functions have the form 

- 1 = -■= ^s'' (68a) 

do 

e'-l = ■■= (68b) 

requiring differentiation of the expressions Eq. (|25p and 
Eq. via the partial differentials {da/dS)d/da and 
{db/dS)d/db, where 

da —a A ^ 

dS " a2 + &2 ' ds" a^ + b^' 

Considering first the fcrmionic case Eq. (j25p . we find 



. f 1 ds eas cot(eas) ebs coihiebs) / , , , ,2 , 2 2 -m^s ,r.n\ 

Ae = — - — coth eos — eos csch eos — cot eas -f eas CSC eas (b + a)s e , (69) 

87r2 7o s3 a2 + &2 V 3 / ' ^ ^ 



for which renormalization only requires subtraction of the logarithmic divergence, since the zero-field constant 
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FIG. 6; The dielectric function for spin-1/2 (left) and spin-0 (right) quantum fluctuations. The value is the correction to the 
Maxwell 1, and as in BI theory, a dominantly magnetic field (below the dark line rising to the right in each plot) gives a 
negative correction suppressing the Maxwell energy-momentum. Contours are provided to display larger values obtained at 
strong electric fields (S < 0). 



is differentiated away. Obtaining the meromorphie ex- 
pansion of the integrand by differentiating the Sitara- 
machandrarao identity Eq. (jA8|) . used above in Eq. (|^ . 
we have the numerically more convenient representation 



s ' [Ka + Kb) ds, 



27r2 a2 + 62 J„ 

kn coth.{kTrb / a) {b / a)csch^ {kirb / a) 



(70) 



(eas)2 — fc27i-2 

/ fcTT coth(fc7ra/6) (a/6)csch2(fc7ra/6) 
Iv((efe)2-t-fc27r2)2 + (e6s)2 + k^n^ 

This form also provides the imaginary part 

f af3l3'^lf k7rp 2/3 ^fc^/?' 
Ime — 777^7 -, — coth— — csch ^, e 



■ 27r /32 -f /3'2 ^ klT 
k—1 



(3' (3' 



(71) 

which is again a reflection the instabihty of strong elec- 
tric fields, as confirmed in figure [7] by its suppression in 
dominantly magnetic fields. The polarization function in 
general field configurations can also be exhibited in the 
quasi-statistical form of Eq. P7|) . but the stronger singu- 
larity in the proper time variable persisting in the 
absence of the m-differentiation makes for simpler 'spec- 
tral' functions 

ln(l-s2+ie) and ln(l + s2) 

for the electric- and magnetic- like integrals Ka and Kt- 
The weak field expansion of the dielectric function can 
be obtained by straightforward differentiation of the ex- 
pansion of the effective action, giving 



90tt SloTT m** 



(72) 



For completeness, we exhibit the dielectric functions for 
magnetic-only 

A,-,B).-?^n,,£<£!±?)fii:^ ,73, 

71" Jo (S + 1) ^ K 



and electric-only 



g(g2 , 2) ^ e-''^^ 

(s2 _ 1)2 ^,27^2 • 

^ ' fc=l 



(74) 



backgrounds, recalling /?' — > nm? /eB and /? -Km? /eE 
in the respective limits. 

For the scalar case Eq. (ps)) . 



e ™- " ds eas csc(eas) ebs csc\i (ebs) , , 
^£=- -rrr^ ^ — T — r5 ^ — - 75) 



. b + a 2 
X { ebs cotli ebs — eas cot eas s 



again renormalized by subtraction of the logarithmic di- 
vergence. The identity used above in Eq. (|64p . provides 
the numerically more convenient representation 



As"" 



4 L roo 

e ah ' 



47r2 a2 + 52 



ds se"" ' {Ka + Kb) (76) 



kn 



if,:=-6^f:(-l)^csch(^)(. 



(b/a) coth(fc7r6/a) 
(eas)2 — k^ir"^ 

k-K 



((e6s)2 -f /!c27r2)2 

(a/6) coth(fc7ra/6)^ 



(e6s)2 + k^TT^ 



12 



io--« 10-"* 10-'^ 10-'^ 10-'- 10-'° 10-* 10-^ 10-* 10--° 10-"* 10-" 10-'* 10-'- 10-'° 10-* 10-^ 10-* 




-10- -1 -10-- -lO"* IQ-* 10-- 1 10- -10- -1 -10-- -10-* IQ-* 10-" 1 lO" 



S [ m*/e- ] S [ m*/e- ] 

FIG. 7: The imaginary part of the dielectric function for spin-1/2 (left) and spin-0 (right) quantum fluctuations. For spin-1/2, 
the dielectric function is consistently positive; however, as indicated on the plot itself in the scalar dielectric function, the 
imaginary part is positive to the left and negative to the right of the overlaid (grey) line at 5 = —V. The vanishing of for 
the anti-self-dual S = —V field configuration is a striking difi^erence from the spinor case, which is suppressed when S > ■ 
Details are suppressed when Ae < 10"'^°. 



displaying the imaginary part 



(77) 



which is positive only for dominantly electric (0 < 7^ < 
—S) fields. The limits of magnetic- and electric-only 
background arc obtained immediately from Eqs. (|73[) and 
(j74p by multiplication with —1/2 and insertion of an al- 
ternating (—1)'^ in the sum. The weakficld expansion is 
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14407r 
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(78) 

Numerical evaluations of the EH corrections to the 
spinor and scalar dielectric functions for general field 
strengths are displayed in figure [S] Ae again refiects the 
unusually square character of the EH integrals, though 
in agreement with Born-Infeld theory, dominantly mag- 
netic fields suppress (Ae < 0) the Maxwell tensor. How- 
ever, the boundary for which this magnetic suppression 
is present differs between fermionic and scalar electrody- 
namics, being approximately S oc in the former case 
and iS cx T-" in the latter. Indeed, the transition to aug- 
mentation (Ae > 0) appears to arise in conjunction with 
the growing imaginary part, as seen in figure [T] 



V. KINEMATICAL EFFECTS OF T 

With numerics providing the magnitude of the induced 
vacuum deformation, we can accurately evaluate physi- 
cal situations in which the modification of the Maxwell 



energy-momentum tensor has observable consequences. 
As may be verified by direct calculation (see also discus- 
sion in Section 8 of |16[ ) , the nonlinear electromagnetism 
(whether BI or QED) behaves in this regard as any other 
nonlinear medium and does not alter the Lorentz force 



(79) 



which is dictated by the necessity of gauge invariance in 
the coupling of EM potentials to charged matter. This 
point is to be contrasted with the modification of particle 
properties. 

Violation of the superposition principle [23, [2^ entails 
an interaction between the background field and the field 
generated by the charged matter. Such an interaction 
could be introduced in the Lorentz force, but it is more 
easily evaluated separately using the weak-field expan- 
sion of the EH effective potential. 

As an example, take a large (r > Ae) charged sphere 
in a strong background magnetic field, which could pro- 
vide a rough model for an a-particle in the atmosphere 
of a highly magnetized neutron star. In the rest frame 
of a non-relativistic charged probe particle, we take the 
background magnetic field as constant B = Bz and the 
electric field as the particle's Coulomb field E = Zer jr^ . 
Integrating the energy of the combined field configura- 
tion, r''^, over the volume with a short distance cutoff 
at the Compton wavelength A, the leading contribution 
is the trace T Wcff : 



"off 



d4^4^(7p2^452) = -^l^(Zei3)2, (80) 



45to4 



45m4 3A 



keeping only the nonlinear-sourced cross terms. The co- 
efficient of the Maxwell energy — also induces cross 
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terms subleading at 0{a^). The cutoff arises since at dis- 
tances shorter than A we must use quantum dynamics to 
describe the probe particle, consideration of which would 
be inconsistent with the classical particle dynamics. 

This interaction energy is positive, independent of the 
sign of the charge, and comparable to the gravitational 
potential of a neutron star with dipolar magnetic field. 
As Ugrav is negative and oc r^^ and the effective (scalar) 
potential goes with oc r^^, 

Meff ^ SnjaZy /eSsurfy i?Lf A.48Af0my ' ^^^^ 

Mgrav 135A V "^e / '''' V / 

converting Newton's constant into the convenient units 
G ~ 1.48 km/solar mass. The A cutoff in Eq. ((80|) can- 
cels against particle mass m, making Eq. ([HT|) indepen- 
dent of both the mass of the particle m and the cutoff 
%. Remarkably, at the surface of a 1.5 M©, 14 km radius 
star with critical surface field -Bsurf = Be, the nonlinear- 
electromagnetic effective potential is 34 times the grav- 
itational potential, resulting in a large repulsive, quasi- 
Lorentz-scalar potential for charged particles entering the 
strong field region. 

For a relativistic particle, the stellar magnetic field 
is Lorentz-transformed requiring consideration of further 
cross-terms. The calculation is simplified by choosing a 
frame, though the effective potential is still determined 
from Eq. (|80|) . 7-factors from the Lorentz transforma- 
tion enhance the effect, but terms linear in the field of 
the particle introduce charge dependence. 



VI. DISCUSSION AND CONCLUSIONS 

The central motivation of this study is the observa- 
tion that externally applied fields in nonlinear electro- 
magnetism have a dark energy-like contribution to their 
energy-momentum tensor. We therefore examined the 
physics giving rise to the trace of the energy-momentum 
tensor as an avenue of insight into the origin of the ob- 
served dark energy in the universe. As in our study is 
generated by quantum-induced nonlinearity of the elec- 
tromagnetic field the physics of dark energy is accessible 
to laboratory experiment probing clcctromagnctism at 
high fields. 

We derived the energy-momentum tensor for general 
nonlinear electromagnetic theories and emphasized the 
form Eq. PH)) . We considered the relationship of the 
trace with the matter condensate and obtained a result, 
Eq. (|36p which amounts to removal of the leading term in 
S in {^Jtp). This reduces the numerical results by a fac- 
tor of about 100 and along with this, the physical dark 
energy effect of the energy-momentum trace is greatly 
reduced. 

Employing the resummation technique introduced 
in [20| , we numerically evaluated the deviations from the 
Maxwell tensor, the dielectric function and the trace for 
Born-Infeld electromagnetism and the Euler-Heisenberg 



effective action with both Fermi and Bose matter fields. 
We believe that these are the first presentations in lit- 
erature of the matter condensate and energy-momentum 
trace arising from the Euler-Heisenberg action at field 
strengths well beyond critical. The dielectric function is 
found in both BI and EH to suppress the Maxwell ten- 
sor in the presence of dominantly magnetic fields. The 
dielectric response of the vacuum thereby enhances the 
observable consequences of the presence of the energy- 
momentum trace. 

The Born-Infeld theory is at first sight an interest- 
ing source of energy-momentum trace. However, a lower 
limit on the limiting BI electric field strength obtained 
30 years ago from the study of precision atomic and 
muonic spectra requires M^/e > 1.710^^V/m, im- 
plying M > 60MeV. Contemporary g ~ 2 experimental 
results, if analyzed with the objective to set a limit on the 
BI scale, would very probably push this limit further up. 
As 5 — 2 of the electron in strong fields is in itsel f a p roject 
requiring study of the two-loop Lagrangian [27|, such 
analysis takes us far beyond the scope of this paper, and 
we leave the investigation to the future. The scale of the 
energy-momentum trace of BI type must be rather large 
and the effect at best comparable to the effect of vacuum 
fluctuations expressed by the Euler-Heisenberg effective 
action. Even though the energy-momentum trace is sup- 
pressed by the QED coupling, (a/vr)^ = 6 x 10~^ (see 
section |IV| . the effective scale 40 Trie — 20 MeV implies 
that quantum fluctuations remain dominant compared 
with any other current theoretical framework, given the 
experimental constraints. 

This is consistent with BI being a high-cutoff theory 
arising from more fundamental matter properties. Hence, 
we devoted the largest part of this report to expand- 
ing our prior study of electron fluctuations in the vac- 
uum, the energy-momentum trace has previously been 
discussed as a signal of vacuum deformation Q . The en- 
ergy density in the trace is then interpreted as the shift 
of the vacuum energy induced by the applied field, and 
in concordance with this interpretation, the trace is pos- 
itive definite when extracted correctly from the effective 
action. The smooth BI T highlights the extraordinary 
form of the quantum-induced trace. Prior considera- 
tions of the analytic structure of the Euler-Heisenberg 
effective action had not made apparent the near singular 
boundaries present in the condensate and the trace. The 
present evaluations suggest further investigations into the 
strong field vacuum phase structure. 

We provided further the first complete numerical cal- 
culations of the electron-positron condensates of spinor 
and scalar QED for arbitrary fields and have addressed 
the contradictory claims in literature relating to claims 
that T and —m{il)'ip) arc equal, and we found in our non- 
perturbative study a clear difference originating in the 
scale dependence of charge renormalization. We empha- 
size that our evaluation of the condensate and of the 
energy-momentum trace are completely independent of 
renormalization procedure. 
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Being a vacuum phenomenon, the QED T studied here 
is a volumetric, extrinsic property, so to obtain observ- 
able effects in the laboratory, one could imagine creat- 
ing a macroscopic region of ultra-strong electromagnetic 
field. On the other hand, as the energy-momentum trace 
of QED is induced by external fields, the effects exist any- 
where and everywhere an electromagnetic field is present. 
For instance, one computes that to generate the observed 
dark energy density 



A 

4^ 



(2.325meV)'' ~ 6.09 x lO'^^^J/m 



(82) 



corresponds to the trace induced by a magnetic field of 
108 T spanning the universe, clearly an unphysical condi- 
tion. The energy-momentum trace induced by nonlinear 
electromagnetism indicates the possibility of observing 
anti-gravity-like effects in the laboratory with high in- 
tensity laser experiments. 

Although many previous studies of vacuum energy and 
issues arising from conformal symmetry breaking have fo- 
cused on QCD, we did not address here how the structure 
of the QCD vacuum, which is very strongly deformed by 
glue and quark fluctuations, relates to the trace Tand 
responds to an applied electromagnetic field. Some dis- 
cussion of this question, including its relation to dark en- 
ergy, has been already offered Combining the quan- 
tum vacuum with general relativity remains today a very 
delicate question however 3. Our report establishes an 
important connection, tying already-recognized quantum 
vacuum effects to dark energy. 

To summarize, we have studied the energy-momentum 
tensor of nonlinear electrodynamics emphasizing an 
explicit relationship of the dark energy-like trace of the 
energy-momentum tensor Eq. (|10b[) to the nonlinearity 
of the theory. In the consideration of electrodynamics as 
a quantum gauge theory, the connection provided a new 
derivation of a non-perturbative identity between the 
energy-momentum trace and the gauge and matter con- 
densates, Eq. ((36)) . The Euler-Heisenberg effective action 
provided a natural example for the numerical evaluation 
of the condensate and energy-momentum trace, the 
results of which are displayed for both fermionic and 
scalar fields in figures [3] and [5] Finally, we briefiy 
explored the implications of an energy-momentum trace 
for charged-particle kinematics. 



APPENDIX A: IMPROVING CONVERGENCE 
OF EULER-HEISENBERG INTEGRALS 

In this appendix, we display the steps in the trans- 
formation of the proper time integrals Eqs. (j^S]) and ([^5]) 
into the more rapidly convergent representations used for 
numerics in the text. 

The effective action for electrodynamics displays non- 
analyticities that, generating an imaginary part of the 
action, are associated with the instability of the vacuum. 
However, our method of resumming the poles is very use- 
ful for improving the overall convergence of integrals of 
the Euler-Heisenberg form, and we start with the case of 
only a magnetic field being present, for which 

f°° ds 

-mM) = — ^ -^{eBscoiheBs-\) (Al) 
Att^ Jo s 

is analytic on the real axis. We use the (subtracted) 
mcromorphic expansions. 



; coth X — \ 



2x 



1 



2 V 

^ -I- /C^TT^ 
fc=l 



oc -. 



Inserting Eq. (|A2|) in Eq. (jAip we obtain 
■n?{eBf 



27r2 



dsY^ 

k=l 



{eBsf + {k-Kf 



(A2) 
.(A3) 

(A4) 



All terms are individually absolutely convergent, so we 
reorder the sum and integral following the procedure 
in (23. After rescaling s skw/eB, the fc-sum is evalu- 
ated in closed form and we obtain 



m{Tjj'ip) 



2^2/3' 



ln(l - e-/^'^) 



ds. 



(A5) 



which is Eq. (|4T|l . 

Eq. (|A3p allows us to remove the quadratic term in 
Eq. (j40p . and, rescaling and resumming, we find Eq. (|42|1 . 
A further integration by parts results in 



27r2 



s — Arctan s 
e^'" - 1 



ds 



(A6) 
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which is the a ^ limit of Eq. p?)) . but we retain the 
form Eq. (|^^ for clarity in the associated discussion of 
signs. 

Turning now to B = 0, i.e. electric field only with 
a ^ \E\, we see in the meromorphic expansion 



X cot X — 1 — 2x 



2 1 

^-^ x"^ — A;27]-2 
k=l 



the singularities that indicate the instability of the sys- 
tem to produce real pairs. We assign to the mass a small 
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imaginary component —^ w? + so that 



m{'ijj'il}) 



27r2 



(A7) 

whence resummation produces Eq. (|43p . Removing the 
leading term in meromorphic expansion for the case of 
the electric field by use of 



X cot X — 1 



k=l 



1 1 



we obtain Eq. 

For general fields, the proper time integrals are rewrit- 
ten using the Sitaramachandrarao identity (Eq. (6) in 



x^ — r, 1 coth(A;7rw/x) 
coth a; cot y = IH '2.x yy , „ „ 



k=l 



kn x'^ + /c^TT^ 



k=l 



k-K — fc^TT^ 



with the result 



(A9a) 



4 3. Z""!, 2 -.V- COth 
•^0 fc=l 



cot h(fc7rVa) 

2_(e»,)2) (A9b) 

coth(A:7ra/6) 



fc7r(fc27r2 4- (e6s)2) 



,(A9c) 



Rescaling converts these expressions to those found in 
Eq. (061). 

The quasi-statistical representation of la in Eq. ((17)) is 
derived by exchanging the sum and the integral in order 
to integrate by parts, 



la^-b ds{s+- In 



1 , /I — X + ie 



2 \l + x - ie 



We break up the coth function: 



coth — 



fe=i 



Since exp(— /s7r6/a) < 1, the n-sum is absolutely conver- 
gent, though slowly when b/a ^ 1. We can exchange the 
order of summation and do the fc-sum: 



1, fl-x + ie\^^ e-P^" 



n,fT 







(A13) 

with f3 = TTm? /ea and m?T-i„ the inner expression in 
Eq. (|A12|) . thus obtaining Eq. (|47)) . 

For the scalar case, Eq. (pS)) . we require identities par- 
alleling those used in the spinor integrations: 



: csch ; 



oc 

1= --x^-2x^y- 



1 (-1)'= 



Ml4a) 



2X-Y.- 



{-If 



fc=i 



^2 + fc27,2 



(A14b) 



X CSC X 



oo 

1- k + 



1 (-1)^ 



2x^^. 



;(Al5a) 



(A15b) 



fc=i 



The representation for general fields uses an identity 
closely related to Eq. (gSl) (see Eq. (20) of [111): 



X — y n \ — ^ 

xy csch X CSC y = l-\ 2x yy 



(— 1)*^ csch(fc7ry/a;) 



k=l 



kn x^ + k^n^ 



'2xy3^ 



(— l)*^ csch{k7Tx/y) 



fe=i 



klT y2 _ /j2^2 



■(A16) 



For the further resummation resulting in the statistical 
representation, the csch function is expanded analogously 



csch X = —e = e ^ > e ^""^ 



1 -e- 



-2x 



n=0 



(A17) 



cothx = > 

^ 1 - 



(AlO) 



(All) 



but the absence of cosh in the numerator means no sum 
over a is introduced. The remaining procedure is the 
same. 



a=±l 



Expanding the denominator as a power series, the sum 
in Eq. (jAlOp becomes 

^ cxp — (m^s -h (2n-H 1 - cr)e6) j . (A12) 
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